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ABSTRACT 

We examine the problem tidally-induced mass loss from collisionless systems such as 
dark matter haloes. We develop a model for tidal mass loss, based upon the phase 
space distribution of particles, which accounts for how both tidal and Coriolis torques 
perturb the angular momentum of each particle in the system. This allows us to study 
how both the density profile and velocity anisotropy affect the degree of mass loss — 
we present basic results from such a study. Our model predicts that mass loss is a 
continuous process even in a static tidal field, a consequence of the fact that mass loss 
weakens the potential of the system making it easier for further mass loss to occur. We 
compare the predictions of our model with N-body simulations of idealized systems in 
order to check its validity. We find reasonable agreement with the N-body simulations 
except for in the case of very strong tidal fields, where our results suggest that a 
higher-order perturbation analysis may be required. The continuous tidally-induced 
mass loss predicted by our model can lead to substantial reduction in satellite mass in 
cases where the traditional treatment predicts no mass loss. As such, our results may 
have important consequences for the orbits and survival of low mass satellites in dark 
matter haloes. 

Key words: gravitation; methods: N-body simulations; galaxies: haloes; galaxies: 
interactions; dark matter 



1 INTRODUCTION 

It is well known that our Galaxy is surrounded by a family 
of satellite systems. These include the several known dwarf 
spheroidal galaxies as well as the Large and Small Magel- 
lanic Clouds and numerous globular clusters (Mateo 1998) . 
The gravitational tidal influence of our Galaxy is currently 
causing the Sagittarius galaxy and globular cluster Palomar 
51 to be torn apart (Ibata et al. 1994; Odenkirhen et al. 
2001). Tidally induced mass loss may also be important for 
explaining the amount and distribution of intracluster light 
(Calcaneo-Roldan et al. 2000; Zwicky 1951). An improved 
understanding of the way in which tidal mass loss occurs is 
therefore of great importance for understanding the evolu- 
tion of galaxies and clusters. 

Furthermore it is known that, although N-body simula- 
tions of cold dark matter (CDM) successfully reproduce the 
overall structure of the Universe, they show potential dis- 
crepancies with observations at galactic scales. The discrep- 
ancy between the number of low mass sub-haloes predicted 
by CDM and the observed number of satellite galaxies is on 
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the order of a factor 5-10 more (Moore et al. 1999; Tasit- 
siomi 2003). This is the so-called satellite catastrophe for 
hierarchical galaxy models. 

Tidal interactions between galaxies may help to resolve 
this apparent discrepancy (Bullock et al. 2000; Benson et 
al. 2002). Sub-haloes can fail to be "lit-up" by a galaxy be- 
cause their star formation has been inhibited by early and 
significant mass loss. There is also the possibility, in extreme 
cases, for the satellite to be totally destroyed due to the grav- 
itational influence of a larger dark matter halo or massive 
galaxy. 

Finally, observations show that merging is a key driver 
in the evolution of galaxy morphologies. N-body simulations 
of the assembly of galaxies in a hierarchical Universe show 
that, as dark matter halos coalesce, the embedded galaxies 
merge on a time-scale that is consistent with dynamical fric- 
tion estimates based on their total mass (Navarro, Frenk & 
White 1995). Accurate treatment of tidal interactions is im- 
portant if we wish to explore this phenomena since dynam- 
ical frictional timescales depend upon the remaining bound 
mass of each satellite. 

Recently, analytic models of satellite orbits and merging 
have been developed and employed to study the properties 
of the sub-halo populations of CDM haloes (Taylor & Babul 
2001; Benson et al. 2002; Benson et al. 2005; Penarrubia 
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et al. 2004; Just & Penarrubia 2004), and to make quan- 
titative predictions for the distributions of such sub-haloes. 
These models typically employ what we will refer to as the 
"classic" treatment of tidal mass loss — namely computing a 
radius in the satellite at which internal and tidal forces bal- 
ance and removing all mass beyond that radius — together 
with simple estimates of the time over which this mass is 
lost. This "classic" method of determining tidal mass loss is 
clearly oversimplified — as particles orbit within the satellite 
halo they may be stripped away by tidal forces as they ap- 
proach apocentre, even though they spend most of their time 
at smaller radii. In this work, therefore, we re-examine the 
process of tidal mass loss and present a more realistic calcu- 
lation based upon the orbits of individual particles within a 
satellite. 

The layout of this paper is as follows. Section [2] briefly 
reviews the "classic" calculation of tidal radius, before de- 
scribing our improved model and presenting some basic re- 
sults from it. Section [3] compares the results of our model 
with N-body simulations. Finally, in [|4] we discuss our re- 
sults and present our conclusions. 



2 ANALYTIC CALCULATION OF 
POST-TIDAL DENSITY PROFILE 

We begin by reviewing the classic calculation of tidal mass 
loss in i\2.1\ before presenting our improved calculation in 

3231 



2.1 Classical Calculation 

In order to determine the tidal radius, we consider a small 
mass orbiting at radius r from the centre of a satellite and let 
the mass of satellite enclosed inside this orbit be M s (< r). 
The satellite itself orbits within a host potential at a radius 
r' . There are two forces on this mass: the gravitational pull 
towards the centre of the satellite and the gravitational tidal 
force from the host potential. 

The gravitational acceleration on the mass towards the 
satellite centre is given by 



GM S y s (x) 



(1) 



where x = r/r s , y s (x) = M s (< r)/M B , and M s and r s are a 
characteristic mass and radius respectively in the satellite. 
The gravitational tidal acceleration on the mass due to the 
host potential can be expressed as: 



GM S d y'(x') 
r 2 dx' x' 2 



(2) 



where x' = r'/r s , y'(x') = M'(< r')/M s , and M'(< r') is 
the mass of the host system within radius r'. Note that we 
have approximated the tidal field as being linear across the 
satellite — the usual approximation made in considering tidal 
fields. The tidal radius is reached when the gravitational 
acceleration on the mass towards the satellite centre equals 
the tidal acceleration on the mass from the host potential, 
i.e. where a s — a t or 



. d y'jx') 
dx' x' 2 



(3) 
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Figure 1. The tidal radius in the classical model plotted as func- 
tion of the tidal field strength, /tidal i f° r the satellite galaxy 
we assume an NFW halo profile with concentration parameter 
C NFW = 11- 



The above equation can be re-expressed in terms of the 
tidal field /tidal, defined to be the gradient of the tidal ac- 
celeration due to the host potential 



/tidal — 

where 

Ctidal = 

giving 



d y'(x') 



dx' 



= Ctidal" 



2/V) 



d y'(x') 



y'{x') dx' x K 



/tidal ■ 



(4) 



(5) 



(6) 



The default values we choose are typical of a massive satellite 
in a Milky Way-like system (M s = 1O 1O M ,M' = 1O 12 M 
and a distance between the centre of mass of the satellite 
and that of the host galaxy of R = lOOkpc, resulting in 
/tidal = 5.8 x 10 -3 — we will explore variations of a factor of 
10 about this default /tidal) - The tidal radius can be found by 
solving the above for x. Figure [T] shows the tidal radius as a 
function of /tidal for a satellite described by an NFW density 
profile (Navarro, Frenk & White 1995) with a concentration 
parameter (the ratio of the virial radius to the characteristic 
scale radius) of cnfw = 11- (For the characteristic radius 
and mass we have chosen r s to be the scale radius of the 
NFW profile and have set M s to equal M200 — the virial mass 
of the halo.) 

2.2 Phase Space Calculation 

As mentioned in the Introduction, the classic model of tidal 
mass loss is clearly highly simplified. It assumes a sharp cut- 
off in the satellite density profile at the tidal radius, and is 
insensitive to the phase space distribution of the satellite 
(i.e. satellites in which orbits are purely tangential would 
have the same tidal radius as those with purely radial or- 
bits, providing that the overall density profile was the same). 
Furthermore the classic model does not take into account the 
orbit of stars, we know that prograde orbits are more eas- 
ily stripped than radial orbits; while radial orbits are more 
easily stripped than retrograde orbits (Read et al. 2006). 
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Figure 2. The geometry of the satellite and particle orbit used 
in our calculation of mass loss. The upper diagram coincides with 
the orbital plane. The particle is identified at position r at which 
it has velocity v. The angle between r and v is <j>. The orbit of the 
particle has pericentre and apocentre r per ; and r apo respectively. 
At apocentre, the particle velocity is v apo . The orbit is assumed 
to precess within the satellite potential (as schematically indi- 
cated by the dashed ellipses) and eventually sweeps out a region 
indicated by the dotted circle. The diameter of this circle is Ar. 
The lower diagram coincides with the plane containing the line 
which connects the satellite centre to the centre of the host po- 
tential, and the normal to the orbital plane. The angle between 
these two lines is £. 



In order to understand the effects of a tidal field on 
the particles (i.e. particles of dark matter or star^J) within 
the satellite with mass M sat , we assume a satellite which 
experiences a constant tidal force (e.g. on a circular orbit 
within a spherical host potential with negligible dynamical 
friction) . Although the assumption of a constant tidal force 
is unrealistic (few satellites will be on orbits that are close to 
circular and, for our default satellite, the dynamical friction 
timescale is of order lOGyr, resulting in significant change 
in orbital radius on a timescale comparable to the mass loss 
timescale from the satellite), we choose to explore this ide- 
alized case here and will explore more realistic examples 
in future work. Particle orbits within the satellite however, 



have arbitrary ellipticity. Figure[2]shows the geometry of the 
particle orbit considered below. 

The energy per unit mass of a particle is 

E(x) = \v 2 + *(*), (7) 

where &(x) is the gravitational potential of the satellite and 
v is the velocity of the particle. We assume an NFW poten- 
tial for the satellite, resulting in a potential: 



$(x) = 



(8) 



where p is the mean density of the halo. 

We consider a rotating frame of reference which rotates 
at the same rate as the orbitting satellite. As the particle 
moves around its orbit it experiences a torque (relative to 
the centre of mass of the satellite) due to tidal (from the 
host potential) and Coriolis forces. The angular momentum 
of the particle therefore varies around the orbit. The parti- 
cle must always remain bound to the combined satellite plus 
host potential (since those potentials are not time-varying 
the energy of the particle is conserved along its orbit). How- 
ever, we can ask the question: If the host potential were 
instantaneously removed, would the particle be bound to 
the satellite or not? To answer this question we simply need 
to know the angular momentum of the particle relative to 
the satellite centre, as this immediately gives the kinetic 
energy of the particle relative to the satellite centre. We ex- 
pect this criterion to be a reasonable discriminator between 
particles which are essentially confined to the satellite and 
those which can move more freely throughout the host po- 
tential (i.e. those in tidal tails). Experiments with N-body 
simulations (see Sj3| confirm this expectation. 

The calculation of the perturbation to the particle's an- 
gular momentum is given in Appendix[X] where we compute 
the maximum value of the particle's energy (in a frame in- 
stantaneously moving with the satellite), _E max . The total 
mass loss suffered from the satellite galaxy is found by inte- 
grating over phase space 



M< 



A-Kr 3 s x 2 F(x, v)27re 2 sin<j!>sin£ 



o J o J a J o 
x H{E max ) df dcpdvdx . 



(9) 



For this work we assume that the phase space density, 
i ? (x,v), is a separable function, so it can be written as: 



4-7rrf p(x)x 2 f(v, 4>)2irv 2 sm<j>sm£ i 



o J o J qJ o 
xH(E ma . x )d£d(f)dvdx 

where p(x) is the density profile of 
1 if s > and H(s) = if s 
the velocity distribution function 
satellite for which we assume an 
Gaussian distribution and which 

Jo Jo ^-O^ 7 ™ 2 sm <j>d<fidv — 1. 
file can be found in a similar way 
simply not performing the integral 



(10) 



the satellite and H(s) = 
< 0. Here, f(v,4>) is 
of the particles in the 
(anisotropic in general) 
is normalised such that 
The radial density pro- 
as shown in eqn. 1101 by 
over radius. 



Our calculation applies to any collisionless particle and can also 
be trivially extended to composite systems, such as a stellar disk 
within a dark matter halo. 



2.2.1 Calculation of the velocity dispersion 

We find the velocity dispersion of the satellite by solving the 
Jeans equation: 
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Figure 3. The evolution of the mass density distribution of the 
satellite according to the model described in EI2.2I The results arc 
shown after 1, 2, 3, 4 and 5 mass loss iterations, and for different 
values of the tidal field: /tidal = 5.8 X 10~ 4 (left panel), 5.8 X 10~ 3 
(middle panel) and 5.8 X 10 -2 (right panel). The initial — prior 
to stripping — satellite profile (purple line) is given for reference. 
Results are shown for an NFW satellite with concentration pa- 
rameter c = 11 and velocity anisotropy parameter f3 = 0. 




Figure 4. The evolution of the satellite's bound mass distribution 
according to the model described in £12.21 where y(x) = Mt,ound(< 
x)/M B , with Mb oun d(< x) being the remaining bound mass of 
the satellite within radius x. The results are shown after 1, 2, 
3, 4 and 5 mass loss iterations, and for different values of the 
tidal field: / tida i = 5.8 X 10" 4 (left panel), 5.8 X 1CT 3 (middle 
panel) and 5.8 X 10 -2 (right panel). Results are shown for an 
NFW satellite with concentration parameter c = 11 and velocity 
anisotropy parameter /3 = 0. In the middle and right-hand panels 
the purple line corresponds to the predictions from the classical 
approach. This method predicts zero mass loss for the plot in the 
left panel as the tidal radius is greater than the virial radius of 
the satellite. 



Figure 5. The evolution of the satellite's bound mass distribution 
(left panel) and mass density distribution (right panel) according 
to the model described in £|2.2l for a velocity anisotropy parameter 
(3 = 0.6. The purple line on the left panel corresponds to the 
predictions from the classical approach, on the right panel the 
initial — prior to stripping — satellite profile (purple line) is given 
for reference. The results are shown after 1, 2, 3, 4 and 5 mass loss 
iterations. The tidal field strength is /tidal = 5-8 x 10 — 3 . Results 
are shown for an NFW satellite with concentration parameter 
c = 11. 
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Figure 6. The satellite's bound mass distribution (left panel) and 
mass density distribution (right panel) according to the model de- 
scribed in i|2.2l for different values of velocity anisotropy parame- 
ter /3 = 0. (red line), /3 = 0.6 (blue line), and /3 = 0.9 (green line). 
The purple line on the left panel corresponds to the predictions 
from the classical approach, on the right panel the initial — prior 
to stripping — satellite profile (purple line) is given for reference. 
The results are shown after one mass loss iterations. The tidal 
field strength is /tidal = 5.8 X 10 — 3 . Results are shown for an 
NFW satellite with concentration parameter c = 11. 



p(x) dx 



p(x)o\(x) + 2 



(3a, 2 (x) _ d§(x) 
x dx 



(11) 



where a r (x) is the radial velocity dispersion. For the radial 



anisotropy (denned as f3 



i/cr r , where o~$(x) is the 



tangential velocity dispersion) we use a fixed value (i.e. in- 
dependent of radius), ranging from (isotropic motion) to 
(3 = 1 (radial motion). Integration is begun at three times 
the virial radius, which is sufficiently large to give an accu- 
rate determination of the velocity dispersion throughout the 
satellite. 



2.2.2 Calculation of the timescale for mass loss 

We wish to estimate the time over which the mass will be 
lost from the satellite. We expect that any unbound particle 
will be lost in a time of order the dynamical time. Defining 
the dynamical time as 



tdyn(x) = 7T 



(12) 



GM s y(x) ' 

an approximate mass loss rate from the satellite is given by 
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Figure 7. The evolution of the satellite's bound mass distribution 
(left panel) and mass density distribution (right panel) according 
to the model described in i|2,2l for a concentration parameter c = 
15. The purple line on the left panel corresponds to the predictions 
from the classical approach. The results are shown after 1, 2, 3, 
4 and 5 mass loss iterations. The tidal field strength is /tidal = 
5.8 X 10~ 3 . Results are shown for an NFW satellite with velocity 
anisotropy parameter (3 = 0. 
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Figure 8. The satellite's bound mass distribution (left panel) and 
mass density distribution (right panel) according to the model 
described in ij2.2l for different values of concentration parameter 
c = 11 (red line), c = 15 (blue line), and c = 20 (green line). 
The purple line on the left panel corresponds to the predictions 
from the classical approach. Results are shown after one mass 
loss iterations. The tidal field strength is /tidal = 5.8 X 10 -3 . 
Results are shown for an NFW satellite with velocity anisotropy 
parameter (3 = 0. 



4nrsX 2 p losa (x)t d y n (x)dx, 



(13) 



where pi OS s(a;) is the density distribution of mass lost 
through tidal stripping (i.e. the original density profile mi- 
nus that after tidal mass loss). The characteristic timescale 
for mass loss is then 



(Wn> = f T M loBB /M lo 



(14) 



where we have introduced a parameter f T to allow this 
timescale to be adjusted to match numerical results. We 
expect f T ~ 1. 



2.2.3 Continuous mass loss 

Using the model described above, we can compute the mass 
lost by the satellite due to tidal forces, and the resulting den- 
sity profile of the remaining bound material. This new den- 
sity profile will give rise to a new potential <I>(r) , weaker than 
that of the original satellite. As such, if we apply our same 
calculation of mass loss to the remaining bound particles, 
using this new potential, we expect that further mass loss 
will occur. This process can, in principle, be repeated ad in- 
finitum. Assuming that the potential changes on a timescale 
of approximately (ta yn ) we can use our model to compute 
mass loss as a function of time. Essentially, we are breaking 
up the continuous process of mass loss into discrete intervals 
(each of order the dynamical time). We compute the total 
mass loss over each interval, before updating the satellite 
density profile and potential prior to computing mass loss 
for the next interval. We will refer to these intervals as "mass 
loss iterations". 



2.3 Limitations of the Model 

Our model, while a vast improvement over the "classic" cal- 
culation of tidal mass loss, has its own limitations: 

(i) In our model there is the implicit assumption that the 
mass loss takes place instantaneously after a time (idyn)- 
Since mass loss is a continuous process our choice to update 
the satellite's potential only after a time (tdyn) is an approx- 
imation. Our model could be modified to use smaller steps 
in time, for example stepping forward by 10% of the mass 
loss timescale, while removing only 10% of the mass pre- 
dicted to have been lost. This would result in more frequent 
updates in the density and potential thereby improving this 
approximation. 

(ii) In computing the energy gained by each particle, we 
evaluate the tidal perturbation by integrating around the 
orbit which the particle would follow in the absence of any 
tidal field. In reality, the orbit will be stretched by the tidal 
field, which will in turn alter the energy gain. These higher 
order contributions to the energy gain are ignored in this 
work, although in principle they could be comparable to the 
first order energy gain. 

(iii) After each mass loss iteration, we hold the density 
profile of the remaining particles fixed and compute a new 
velocity dispersion profile using eqn. pip . This effectively 
changes the energies and angular momenta of the remaining 
bound particles. A correct calculation of the phase space 
distribution of the particles after each episode of mass loss 
is an important missing ingredient of our model. We defer 
detailed study of this issue to a future paper (see Sj4] for 
further discussion of this point). 



2.4 Basic Results 

In this section will present the basic results of the model 
discused in t|2.2l Unless otherwise stated, we will consider 
an NFW satellite with concentration parameter cnfw = 11, 
isotropic orbits (/3 = 0) in an /tidal = 5.8 x 10~ 3 tidal field. 

Figure 3 shows the evolution of the mass density dis- 
tribution of the satellite. The results are shown after 1 (red 
line), 2 (blue line), 3 (green line), 4 (orange line) and 5 (cyan 
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line) mass loss iterations, and for different values of the tidal 
field: / tida i = 5.8 x fCT 4 (left panel), 5.8 x 1CT 3 (middle 
panel) and 5.8 x 10~ 2 (right panel). The initial — prior to 
stripping — satellite profile (purple line) is given for reference. 
As expected, stronger tidal fields result in greater mass loss 
from the satellite. For the strongest tidal field shown, sig- 
nificant mass loss occurs even at the scale radius (x = f ) of 
the satellite. 

Figure 4 shows the corresponding evolution of satellite's 
bound mass distribution. The results are shown after 1 (red 
line), 2 (blue line), 3 (green line), 4 (orange line) and 5 (cyan 
line) mass loss iterations, and for different values of the tidal 
field: / ti dai = 5.8 x 10" 4 (left panel), 5.8 x f(T 3 (middle 
panel) and 5.8 x 10~ 2 (right panel). In the middle and left 
plots the purple line corresponds to the predictions from the 
classical approach. This method predicts zero mass loss for 
the plot in the left panel as the tidal radius is larger than the 
virial radius of the satellite. Note that our model predicts 
that the satellite continues to lose mass in each successive 
mass loss iteration — there is no evidence that the mass is 
converging to some fixed value. This is despite the fact that 
the tidal field is constant and is due to the weakening of 
the satellite's potential after each mass loss iteration. The 
classic model for tidal mass loss predicts a fixed mass for the 
satellite after tidal stripping. In our model, the satellite falls 
below this mass after two mass loss iterations for (except 
for the weakest tidal field shown, in which case the classic 
model predicts no mass loss). 

In Figs. 5 and 6 we explore the dependence of mass loss 
on the velocity anisotropy parameter (5. More specifically, 
Fig. 5 shows the evolution of the satellite's bound mass dis- 
tribution (left panel) and mass density distribution (right 
panel) for a velocity anisotropy parameter /3 = 0.6. The 
purple line on the left panel corresponds to the predictions 
from the classical approach, on the right panel the initial — 
prior to stripping — satellite profile (purple line) is given for 
reference. The results are shown after 1 (red line), 2 (blue 
line), 3 (green line), 4 (orange line) and 5 (cyan line) mass 
loss iterations. 

Figure 6 shows the evolution of the satellite's bound 
mass distribution (left panel) and mass density distribution 
(right panel) for different values of velocity anisotropy pa- 
rameter f3 = (red line), /3 — 0.6 (blue line), and /? = 0.9 
(green line). The purple line on the left panel corresponds 
to the predictions from the classical approach, on the right 
panel the initial — prior to stripping — satellite profile (purple 
line) is given for reference. The results are shown after one 
mass loss iteration. From Figs. 5 and 6 it can be seen that 
increasingly radial orbits result in less mass loss. 

We note that our current calculation do not accurately 
account for the ellipticity of orbits (see Appendix . As 
such, we must treat this result with some caution. A more 
detailed analysis is deferred to a future paper in which the 
shapes of orbits will be accounted for. 

In Figs. 7 and 8 we explore the dependence of mass 
loss on the concentration parameter cnpw- Figure 7 shows 
the evolution of the satellite's bound mass distribution (left 
panel) and mass density distribution (right panel) for a con- 
centration parameter cnpw = 15. The purple line on the 
left panel corresponds to the predictions from the classical 
approach. The results are shown after 1 (red line), 2 (blue 



line), 3 (green line), 4 (orange line) and 5 (cyan line) mass 
loss iterations. 

Figure 8 shows the evolution of the satellite's bound 
mass distribution (left panel) and mass density distribution 
(right panel) for different values of concentration parameter: 
cnfw = 11 (red line), cnpw = 15 (blue line) and cnfw = 20 
(green line). The purple line on the left panel corresponds to 
the predictions from the classical approach. The results are 
shown after one mass loss iteration. From Figs. 7 and 8 it 
can be seen that more concentrated halos experience greater 
mass loss. Note, however, that we keep /tidal fixed for each 
satellite considered here and use the NFW scale length as 
our unit of length, r s . As a result, halos with higher values 
of cnfw actually have a smaller fraction of their total mass 
within x — 1 and are more extended (i.e. their virial radii lie 
at larger values of x) making them more susceptible to mass 
loss. This does not necessarily reflect the scalings expected 
for cosmological halos for which the concentration and mass 
are correlated (Bullock et al. 2001). 



2.5 Comparison to Previous Work 

Here we will compare our method with the results presented 
in (Hayashi et al. 2003). In particular, they found that 
one can get the satelite's density profile after tidal strip- 
ping starting with an NFW profile and taking into account 
two kind of modifications: a lowering of the central density 
/ and the introduction of a tidal radius rt e of the outer pro- 
file. Hayashi et al. provide fitting curves both for the tidal 
radius r te and the lowering of the central density factor / as 
a function of the satellite bound mass, so we are able to see 
if our results agree or not. 

We choose as a comparison variable the tidal radius rte 
since it is a better defined quantity. Since in our model we do 
not have a sharp cut-off radius we define r te as the radius 
at which the density drops to 50% of the original NFW 
value. Figure 9 shows the tidal radius rt e as a function of 
the satellite bound mass. The solid line corrsponds to the 
fitting curve provided by Hayashi et al., the coloured points 
are results from our model that correspond to tidal fields of 
different strength(red points, /tidal = 5.8 x 10 -2 , cyan points, 
/tidai = 5.8 x 10~ 3 , and purple points, / tidal = 5.8 x 10" 4 ). 
Our results, especially for the case of strong tidal filds are 
in good agreement with Hayashi et al. fitting curve. 



3 COMPARISON WITH N-BODY 
SIMULATIONS 

In order to test the ability of our model to quantitatively 
describe mass loss from dark matter haloes we have run 
several N-body simulations. Simulations were performed us- 
ing Gadget2 (Springel 2005), modified to include a point 
mass at the centre of the coordinate system. To study mass 
loss, an N-body representation of an NFW dark matter halo 
with a concentration parameter of 11 was placed in a circu- 
lar orbit around this point mass. A softening length equal 
to 14% of the scale length was chosen (this is the opti- 
mal softening length suggested by Power et al. 2003). To 
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Figure 9. The tidal radius rte as a function of the satellite bound 
mass. The solid line corrsponds to the fitting curve provided by 
Hayashi et al., the coloured points are results from our model 
that correspond to tidal fields of different strength(red points, 
/tidal = 5.8 X 10 -2 , cyan points, /tidal = 5.8 X 10 -3 , and purple 
points, / tidal = 5.8 X 1CT 4 ). 

ensure accurate evolution^, we set the maximum timestep 
to be 13% of the dynamical time at the scale radius (cor- 
responding to 1% of the dynamical time at the virial ra- 
dius) and set ErrTolIntAcc= 0.025. A force accuracy of 
ErrTolForceAcc= 0.005 was chosen. 

We have performed numerous tests to ensure that our 
simulations are not affected by the choice of numerical pa- 
rameters in Gadget2. Key numerical parameters such as 
the softening length, maximum timestep, integration accu- 
racy and force accuracy were reduced by factors of two and 
the calculations repeated. In no case do we see any signifi- 
cant change in the evolution of the satellite mass or density 
profiles. As such, the mass loss seen in these calculations 
is due entirely to the tidal field from the point mass. Fur- 
thermore, with these parameters, the orbital radius of the 
centre of mass of the satellite varies by only 0.06% through- 
out the duration of our simulation. This corresponds to a 
0.18% variation in /tidal- 

The mass and orbital radius combination of the circular 
orbit were chosen to produce a desired /tidal- The satellite 
was allowed to orbit for a period corresponding to over 85 
dynamical times at the scale radius (or, equivalently, almost 
7 dynamical times at the virial radius). At each timestep 
of the calculation we performed a binding energy analysis 
on all satellite particles which were bound at the end of the 
preceding step to check whether any of these had become 
unbound. (We define "bound" here to mean that the parti- 
cle's energy in a frame coinciding instantaneously with the 
centre of mass of the satellite, and including gravitational 
energy only from other bound particles, is negative.) The 
remaining bound particles were used to compute a density 
profile and total mass for the satellite. 

Satellites were set up using 10 5 particles sampled from 
the phase space density described in i|2.2l and with particles 
out to three times the virial radius. 

To test for numerical convergence, we began by running 

2 The reader is referred to the GADGET2 user guide for definitions 
of numerical parameters. 



our calculations using smaller numbers of particles. The par- 
ticle number was gradually increased until converged results 
were obtained. Using 3 X 10 4 particles, for example, we find 
that the remaining bound mass of the satellite differs from 
that in the 10 5 particle simulation by less than 2% at all 
times. As such, 10 5 particles are enough to ensure that par- 
ticle number does not significantly affect our results. 

Figure [10] shows the mass of our satellite in an /tidal = 
5.8 x 10 -3 tidal field as a function of time. The solid black 
line indicates the remaining bound mass of the satellite 
within the virial radius. For reference, the dotted black line 
indicates the total bound mass of the satellite (i.e. includ- 
ing particles beyond the virial radius) — unsurprisingly mass 
beyond the virial radius is quickly stripped away. At several 
times during our simulation we have extracted the parti- 
cles which remain bound at that time and evolved them in 
isolation (i.e. without the presence of unbound particles and 
without any tidal field applied) . Results from three such cal- 
culations are shown as solid grey lines in Fig. 1101 With the 
tidal field removed, the mass of the satellite quickly reaches 
a constant valu^E providing further evidence that numerical 
effects do not contribute to the mass loss from this satellite, 
which must instead be entirely due to the applied tidal field. 
The fact that the satellite mass remains close to that at the 
time at which the host potential was removed validates our 
criterion for deciding whether or not a particle should be 
considered to be bound to the satellite or not. 

Inspection of the N-body mass loss curve indicates that 
there are two distinct regimes of mass loss: an initial rapid 
regime in which mass declines almost exponentially with 
time (illustrated by the thin dashed line) followed by a slower 
regime in which mass declines as a power-law in time (illus- 
trated by the thin dotted line). 

The results of analytical calculations of mass loss are 
also shown in Fig. 1101 The horizontal dashed line indicates 
the mass of the satellite after tidal mass loss according to 
the classical model. The cross on this line indicates the mean 
dynamical time for mass loss in this model for / Tdv „ = 1 
(note that we could always force this cross to lie on the N- 
body line in the classical model by suitable choice of /-r dyn )- 
The limitations of the classical model are clearly seen — the 
satellite loses mass continuously throughout the calculation, 
while the classical model predicts a fixed mass. 

Circles indicate the prediction from the mass loss model 
described in this work. We show results for 17 mass loss it- 
erations and chose /r dyn = 0.35. As expected, /r dyn is of 
order unity. Note that our model is able to match the rate 
of mass loss quite well for about five iterations, before be- 
ginning to significantly overestimate the rate of mass loss. 
This discrepancy will be discussed further in 

Figure [TJJ shows mass loss as a function of time for 
the same satellite in tidal fields of three different strengths. 
Coloured lines indicate the bound mass of the N-body satel- 
lite, while coloured points indicate the results from the 
model described in this work. Dashed lines with crosses in- 



3 Note that the mass does decline slightly even after the tidal 
field is switched off. This transient mass loss is due to the fact 
that the satellite must relax to a new configuration after the tidal 
potential is removed, leading to some particles gaining energy and 
becoming unbound. 
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Figure 10. The mass, y(t) (defined as the remaining bound mass 
divided by the initial bound mass of the satellite), of a c = 11 
satellite in an 

/tidal = 5.8 X 10~ 3 tidal field as a function of time 
(in units of the dynamical time at the satellite scale radius). The 
solid black line shows the remaining bound mass of the satellite 
within the virial radius as determined by N-body simulation. The 
dotted black line indicates the total bound mass of the satellite 
(i.e. including particles beyond the virial radius). Grey lines in- 
dicate how the mass of the satellite evolves when the tidal field 
is removed at various points in time. The horizontal dashed line 
indicates the mass within the classic tidal radius, while the cross 
on this line indicates the mass loss timescale (for / Tdyn = 1) for 
the classic calculation. Points indicate the mass predicted by the 
model described in this work after 1-17 iterations. Thin dashed 
and dotted lines are shown to illustrate the two distinct regimes of 
mass loss seen in the N-body simulation — an initial rapid regime 
in which mass declines almost exponentially with time (dashed 
line) followed by a slower regime in which mass declines as a 
power-law in time (dotted line). 



dicate the results of the classical model as in Fig. [10] To 
obtain a good fit to the mass loss rates we are forced to 
adopt different values of f T depending on the strength of 
the tidal field. The results show have f T = 1.4, 0.35 and 
0.15 for tidal fields of /t tda i = 5.8 x 10~ 4 , 5.8 x 10" 3 and 
5.8 x 10 -2 respectively. Once again we see that our model 
describes the mass loss quite well for several iterations be- 
fore overestimating the mass loss rate at late times. (This 
overestimation is not seen for the /tidal = 5.8 x 10~ 2 calcula- 
tion as the satellite loses its mass so rapidly that the second, 
slow mass loss regime is never seen.) 

In addition to the total bound mass of the satellite, 
our model is able to predict the radial density profile of 
that mass. In Fig. [12] we show as coloured lines the den- 
sity profile (normalized to the original density profile) of the 
same satellite in three different tidal fields after 1-5 iter- 
ations. Coloured circles show the density profile from the 
N-body simulation at the corresponding times. Vertical dot- 
ted lines indicate the softening length in our calculations. 
As expected, in the N-body simulations the density at radii 
less than a few times the softening length drops quickly. The 



Figure 11. The mass of a c = 11 satellite in tidal fields of varying 
strengths as a function of time (in units of the dynamical time at 
the satellite scale radius). Solid lines show the remaining bound 
mass of the satellite within the virial radius as determined by N- 
body simulation. The horizontal dashed lines indicate the mass 
within the classic tidal radius. Points indicate the mass predicted 
by the model described in this work after 1-5 iterations. 



vertical dashed line indicates the tidal radius in the classic 
model (this line cannot be seen in the upper panel as the 
classical model predicts a tidal radius beyond the virial ra- 
dius for this tidal field strength). 

For the weakest tidal field (upper panel of Fig. I12J1 the 
analytic model predicts the radial density profile seen in 
the N-body simulation reasonably accurately for the first it- 
eration (red line). After this, our model underpredicts the 
mass loss rate and, as such, overpredicts the resulting den- 
sity profile. As the tidal field is increased the model fairs 
worse. For the strongest field considered we can see that the 
N-body satellite quickly develops a constant density core, 
the density of which decreases with time. Our model does 
not reproduce this behaviour. This may reflect the fact that 
our model applies only first order perturbations to the parti- 
cle orbits. As /tidal is increased the orbits of particles which 
remain bound to the satellite become ever more perturbed 
due to the tidal force. A measure of this perturbation can be 
constructed by averaging —AE/E (see eqn, IA42l for the def- 
inition of AE) over all particles which remain bound, which 
measures the fractional change in particle orbital energies. 
This quantity, after one iteration, is 0.08, 0.21 and 0.30 for 
/tidal = 5.8 x 10 -4 , 5.8 x 10~ 3 and 5.8 x 10" 2 respectively. 
Thus, for the strongest field considered this factor is be- 
coming quite large. Furthermore, for an NFW potential, the 
gravitational potential varies only very slowly with radius 
at radii less than the scale radius. Thus, a relatively small 
perturbation in the energies of particles at these radii can 
lead to a large perturbation in their apocentric distance. As 
a result, the second order correction to the work done by 
the tidal field can be large. 
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Figure 12. The density profile of the satellite galaxy is shown 
relative to the original profile and as a function of radius (in units 
of the satellite scale radius). The three panels show results for 
three different values of /tidal : 5.8x 10 -4 (upper panel), 5.8x 10 — 3 
(centre panel), 5.8 X 10 -2 (lower panel). Lines show predictions 
from the calculations described in this work, while circles indicate 
results from the N-body simulations. Results are shown after 0—5 
iterations as indicated by colour (black, red, blue, green, magenta, 
cyan). Vertical dotted black lines indicate the softening length in 
the N-body simulations. Vertical dashed black lines indicate the 
classic tidal radius. 



4 DISCUSSION AND CONCLUSIONS 

We have presented a model of tidal mass loss from collision- 
less systems with arbitrary phase space distributions. Our 
calculation has many advantages over the classic model of 
tidal mass loss (in which the density profile is truncated 
beyond a tidal radius determined through force balance ar- 
guments). In particular, our calculation takes into account 
the orbital structure of the system, permitting the effects of 
anisotropic orbits on the degree of mass loss to be investi- 
gated. Furthermore, we are able to estimate the density pro- 
file of the material remaining after mass loss has occurred. 

A key prediction from this model is that mass loss will 
be continuous even in a static field — the bound mass of the 
system shows no sign of converging to a fixed value. This 
behaviour is also seen in the N-body simulations which we 
have carried out. In particular, the N-body simulations show 
evidence for two distinct regimes of mass loss: an initial rapid 
phase in which the mass declines exponentially with time 
followed by a slower phase during which the mass declines 
as a power-law in time. 

The rate of mass loss predicted by our model is com- 
parable to that seen in the N-body simulations during the 
initial, rapid mass loss phase, although we find that the pa- 
rameter f T (which scales the rate of mass loss in our model) 
must vary with the tidal field strength in order to match 
N-body results. We intend to explore and characterize this 
systematic variation of f T in a future paper. This continuous 
mass loss is a consequence of the fact that mass loss weakens 
the gravitational potential of the system thereby allowing 
particles to become unbound that could not have escaped 
from the original potential. At larger times, while the N- 
body simulations transition to a slower mass loss regime our 
model continues to show a rapid, exponential mass loss. The 
origin of the slow mass loss regime in the N-body simula- 
tions and the failure of our model to reproduce this is not 
understood at present, but is the focus of ongoing study. 
Nevertheless, our model is able to describe the rate of mass 
loss quite well over almost 10 satellite dynamical times. 

The model presented in this work predicts small 
amounts of mass loss even when the tidal field is sufficiently 
weak that the classic model predicts zero mass loss (due 
to the tidal radius lying beyond the outer radius of the sys- 
tem) . As this mass loss occurs continuously it can eventually 
lead to a large reduction in the mass of the system (e.g. in 
the left-hand panel of Fig. [3] with /tidal = 5.8 x 1(F 4 , the 
mass has been reduced to almost 60% of its original value 
after five mass loss iterations). This could have important 
consequences for the distribution and survival of low-mass 
satellite systems in dark matter haloes. 

In Fig. Q2] we see that, for the strongest tidal field, our 
model fails to predict the reduction in density in the centre 
of the satellite. This may explain the rather low value of f T 
we require to fit the mass loss rate in this case. If our model 
correctly predicted the mass loss from the centre of the halo 
it would naturally predict a lower (r). As noted in £[3] this 
may be due to the fact that our estimate of the tidal torque 
exerted on each particle is correct only to first order — an 
improved calculation accounting for the distortion of particle 
orbits by tidal torques should result in greater mass loss, and 
is expected to be particularly important at small radii. We 
defer study of this issue to a future paper. 
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We have demonstrated that our model is able to con- 
sider the effects of orbital anisotropy and profile shape on the 
degree of mass loss. A full study of the effects of anisotropy 
must await a full treatment of elliptical orbits, to be pre- 
sented in a future paper. 

The major shortcoming of the present model is the lack 
of a reasonable method to compute the phase space density 
of particles after each mass loss iteration. Currently, we sim- 
ply assume that the density profile of the remaining bound 
material remains fixed (i.e. the radial distribution of bound 
particles after mass loss is identical to the radial distribu- 
tion of the same particles prior to mass loss). We then simply 
re-solve Jeans equation to find a velocity dispersion which 
results in an equilibrium system built from these particles. 

A possible improvement upon this may be to treat 
the potential of the satellite as varying only slowly with 
time. Changes in particle orbits can then be considered in 
terms of adiabatic invariants and a new phase space den- 
sity constructed which explicitly conserves these invariants 
for bound particles after mass loss has occurred. It should 
be noted, however, that mass loss is assumed to occur on a 
time comparable to the local dynamical time in the satellite 
(an assumption confirmed by the N-body simulations car- 
ried out here). As such, the changes in the potential occur 
on a timescale comparable to the orbital times of particles, 
making the adiabatic approximation a poor one. The success 
of such an approach therefore remains to be determined. 

In a future paper, we will develop this model further. 
Along with the improvements mentioned above, we intend 
to examine more carefully the dependence of the energy 
gain through tidal torques on orbital eccentricity and the 
effects of using smaller time intervals for each mass loss it- 
eratiorQ. We will also explore the application of our model 
to two-component (i.e. dark matter + stellar) systems and 
present fitting formula for mass loss as a function of tidal 
field strength and time. A further important consideration 
will be to apply our model to situations in which the tidal 
field varies with time. 

The model of tidal mass loss presented here represents 
part of an ongoing effort to improve the accuracy and pre- 
dictive power of analytic models of satellite orbits and evo- 
lution. Detailed modelling of this sort is required in order to 
understand the substructure of dark matter halos and the 
evolution of the galaxies within them. 
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APPENDIX A: ANGULAR MOMENTUM 
PERTURBATIONS 

This Appendix details our calculation of the effects of tidal 
forces on particles orbitting in a satellite based upon per- 
turbations to a particle's angular momentum. 



Al Tidal Torques 

We consider a satellite galaxy placed in a circular orbit 
within a spherical host potential. Effects of dynamical fric- 
tion are neglected. We consider a frame which rotates at the 
same rate as the satellite orbits. 

The gravitational force (per unit mass) at a position r' 
from the centre of the host potential is 



GM'(d) . 
Fg( d ) = j 2 — d > 



(Al) 



where M'(d) is the mass of the host contained within radius 
d and a hat indicates a unit vector. 

In addition, in this rotating frame, the particle experi- 
ences a fictitious force of 



fi d — 2f2 A v, 



(A2) 



where $1 is a vector with magnitude equal to the angular 
frequency of the satellite orbit, and normal to the orbital 
plane, v is the velocity of the particle. For the circular orbit 
considered, 



fi = ^GM'(d)/d s . 



(A3) 



The centrifugal force here will act to cancel the mean grav- 
itational force on the satellite since 



d 2 



(A4) 



At any point in the satellite, we are therefore left with a tidal 
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force (i.e. the gravitational force minus this mean gravita- 
tional force) which has a quadmpole form, the Coriolis force 
and an internal force due to the mass of the satellite itself. 
This latter will be ignored since it exerts no torque on the 
particle. 

We consider a particle orbitting within the satellite. Let 
the satellite centre be at r', and the particle be at d. The 
orbital position of the particle relative to the satellite centre 
is r = d — r'. We are interested in the angular momentum 
of the particle around the satellite centre and so wish to 
compute the torque exerted on the particle relative to this 
centre: 



T = rA[F g (d) + F S2 
= r A 



(A5) 



Considering just the tidal part of the above, the magnitude 
of the tidal torque is 



T tid 



GM'(d) . GM'jr') . ' 

d 2 Sm ff ^— Sm9 



(A6) 



where 9 is the angl^fl between r and r' and fi is the angle 
between r and d. Defining 

„'3 



Ctidal 



GM'(r') 



GM'(r') dr' \ 
(e.g. Ctidal = —2 for a point mass), we can write 



GM'jd) _ GM'(r') 
J2 ~ ~2 



1 + c tid 



d- 



(A7) 



(A8) 



which is the usual linear approximation for the tidal field, 
valid providing r <C r' '. Therefore, 



idal 



GM'(r' 



d — r' 

1 + Ctidal ; — ) sin P - sin f 



If a is the angle between r' and d, then 

sin a sin 8 

r d 
and P = 7r — 9 — a. We also have 

d 2 = r'' 2 + r 2 — 2r'r cos 9. 

Combining these results, we find 

GM'(r') 



(A9) 

(A10) 
(All) 



1 + Ctid 



(0 + sin" 1 j^sintfj) - sin 6 



(A12) 



Substituting for d and expanding as a series in r/r' we find 



Ttidai = r 



GM'(r') 



(1 - ctidal) cos 9 sin 9 — 
r 



1 2 r 

+ -(3cos 9 - 1)(1 - Ctidal) sin 0— 
I r 



(A13) 



The first term in this expression shows the quadrupole na- 
ture of the tidal torque. Note that the tidal torque depends 

5 Note that r A r' and r A d point in the same direction. 



on Ctidal as expected, but there is an additional contribution 
(i.e. we have 1 — Ctidal rather than just —Ctidal) which arises 
from the fact that, even if there is no gradient in the gravi- 
tational force as a function of distance from the host centre, 
there is still a difference in the vector forces acting at the 
satellite centre and at the position of the particle, which acts 
as a tidal field. We will ignore higher order terms from now 
on. 

We now wish to find the change in the angular mo- 
mentum of the particle as it moves around its orbit in the 
satellite galaxy. Note that the tidal torque always acts nor- 
mal to the plane containing r' and r. The change in angular 
momentum around the orbit is given simply by 



Aj 



Ttidaldt. 



(A14) 



Writing dt — &x/~X. where \ is an angle measured around 
the orbit, and using the fact that \ ~ jo/ r (x) 2 j where 
jo = rv sin <f) is the unperturbed angular momentum of the 
particle, we have 



Aj 



1 tidal : OX- 
JO 



(A15) 



For an orbit in a plane whose normal makes an angle £ with 
r', we have 



cos 9 = cos x sin 



(A16) 



where we have chosen \ to coincide with 9 for the case £ = 
7I-/2. Thus 



Aj 



GM'(r') „ 

. ,3 (1 - Ctidal) 



Tcosxsin^l - cos 2 x sin 2 C) 1/2 r 4 (x)dx- (A17) 



To assess the vector change in jo it is useful to consider 
three orthogonal unit vectors defined by f', (f' A jo)/ sin £ 
and [f A (f A jo)]/ sin£. The angular momentum vector jo 
has projections onto these three vectors of 



(A18) 




The tidal torque unit vector has projections 
( 

T = sign(sin6lcos6>) |cos£cos#| 

\ (1 - cos 2 £ cos 2 0) 1/2 

where 



(A19) 



(A20) 



and we have made explicit the quadrupole form of the 
torque. Therefore, the change in the angular momentum is 



-1 


if 


X 


< 





if 


X 


= 


1 


if 


X 


> 



Aj 



GMV) (1 
jor'" 



(1- f'tidai) / r (x)sign(sin0cos( 

2 Ml/2 



x cos x sin £(1 — cos x sm 


/ 1 I cos x cos f sin £| | dx- (A21) 

|(1 - cos 2 ^ cos 2 x sin 2 ^) 1/2 1 
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Even for a circular orbit (i.e. constant r) the integrals are 
not analytically tractable in the general case. For £ = 7r/2, 
we find the change in angular momentum between xi and 
X2 to be 



Aj 



i GM'(r') 

j or /3 



(1 - Ctidal) 



l/2[cos 2 X \ll 



(A22) 



and so, the maximum variation in j is given by (e.g for 
Xl = and X2 = 7r/2) 

4 GM'(r') 



j = jo ± r 



"(I - Ctidal). 



In general we can write 



(A23) 



Aj 



JOfftidal 



r(x) 



sign (sin 6 cos ( 



,1/2 



x cosxsin£(l — cos x sm 


1 | cos x cos £ sin £| 

|(1 - cos 2 ^ cos 2 x sin 2 ^) 1/2 1 



d X . (A24) 



where the dimensionless number gtidai = (GM'(r')r'/jo)(h 
ctidal)- Then, the perturbed angular momentum at x is 



j(x) =jo+Aj(x). 



(A25) 



Considering the kinetic energy (per unit mass) of the particle 
in a frame moving with the satellite centre, 

ij 2 (x) 



2 V{X) 2r 2 (x)' 
we find 

\ v2 (x) = \yl{x) 



1 + 2 Ai_j£ + ^! 

Jo Jo 



and so the gain in energy is 

" Aj( X )-jo , Aj 2 (x) 



j! 



2r 2 ( X ) 

Simplifying to circular orbits 

" Aj(x) - jo , Aj 2 (x) 



AE =-4 



2r 2 



Jl 



(A26) 



(A27) 



(A28) 



(A29) 



Integrating from x = to x = t/2 gives the maximum 
change in kinetic energy (as this maximizes both terms in 
the above). Defining, 

hi(0 = 2Aj(0,^/2)-jo/j ( Vfftidai/(r/r') 4 (A30) 

MO = Aj 2 (0^/2)/j 2 /< ? 2 dai/(r/r') 8 , (A31) 

numerical integration leads to the forms shown in Fig. IA1I 
Writing out the energy expression in detail we find 



AE 



1, , tu1 .GM'(r') 

-/ll (0(1 - Ctidal) 



1, ,+s„ x2 f GM'(r') 

+ 2^(0(1 - Ctidal)' I ^ 



or 



v' 2 / r \ 2 



+ Me)(l- Ct Mair [fj {^J 



(A33) 




Figure Al. The function fci(§) for the cases where Coriolis 
torques are not included is shown by the solid line. Results for 
the case where Coriolis forces are included (and are averaged over 
7) are shown by the dashed line. 



where v' = y GM' (r'JJr' and j' — v'r' . For cosmological 
halos we would expeclu v oc r and so j oc r 2 . As such, since 
r <C r', the second term in the above will be negligible and 



AE = 



T (£) 2 M£)(i-c« d 



(A34) 



For non-circular orbits, the problem becomes more com- 
plicated. The maximum energy changes depends on the 
value of x at which apocentre and pericentre occur, and on 
the shape of the orbit. Experiments using elliptical orbits of 
eccentricity e show that using the semi-major axis in place 
of r and 



(A35) 



is a reasonable approximation (better than 20%) for e < 0.8. 
Application to the more complex orbits found in typical dark 
matter halos is deferred to a future paper. As such, we will 
ignore any dependence of h\ (£) on orbital shape in this work. 

A2 Coriolis Torque 

In the rotating frame, the particle experiences a Coriolis 
force in addition to the tidal force, which can be of compa- 
rable magnitude. This results in a torque 



Tcorioiis = -2r A (O A v) 

= -2[(r-v)n-(r-0)v]. 
For circular orbits r ■ v = 0, so 
Tcorioiis = 2(r • n)v. 



(A36) 
(A37) 



6 At the virial radius, r, of any halo of mass M we have v 2 = 
GM/r. Assuming that cosmological halos at any given redshift 
all have the same mean density within their virial radius, then 
Mar 3 and so v 2 oc r 2 . 
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Note that vector fi is normal to r'. Denning angle 7 as a 
measure of the rotation of the normal to the orbital plane 
around r' with 7 = corresponding to the plane containing 
r' and ft we have 



- Coriolis 



2r Q (cos £ cos 7 cos x — sin 7 sin x) v - ( A38) 



The change in angular momentum due to the Coriolis force 
is then 

Ajcorioiis = / 2—r- (cos£cos7cosx-sin7sinx)vdx-(A39) 
J Jo 

Projecting this onto our basis vectors we find 

Aj = 2^ — ^ v'r' J (cos ^ cos 7 cos x — sin 7 sin x) 

(sin x sin £ \ 
cosx d X . (A40) 
sin x cos £ / 

where we have used v = ja/r and = v'/r'. Alternatively, 

Aj = - gti ^ i j J (cos g cos 7 cos x - sin 7 sin x) 

(sin x sin £ \ 
cosx dx- (A41) 
sin x cos £ / 

Note that, for cosmological halos, we expect character- 
istic properties of halos of different masses to obey v oc r. 
As such, we find gtidai oc (r'/r) 4 . It can then be seen that 
the change in angular momentum due to the tidal force and 
the Coriolis force are of comparable magnitude. 

We can compute the functions hi(£) and /12(C) when 
the Coriolis torque is included, by averaging over angle 7. 
The results for hi (£) is shown as a dashed line in Fig. IA1I 
We find that the inclusion of the Coriolis torque significantly 
alters h\(£) (by up to a factor of around 2). We therefore 
choose to include the Coriolis term in all calculations carried 
out in this work. 

The case of non-circular orbits including Coriolis 
torques becomes significantly more complicated and a full 
treatment is deferred to a future paper. 



A3 Final Result 

The final statement is that the change in energy of the par- 
ticle can be approximated by 

'2 2 

^ E = V ~y(j) hi(0(l-ctid*), (A42) 

where hi(£) is computed with the inclusion of the Coriolis 
torque term as shown in Fig. lAll Thus, the maximum energy 
of a particle relative to the satellite centre is 

S m ax = E + AE. (A43) 

It is this value of -E max which is used in eqn. (|10[) . 



